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■^ ' Abstract 

\^ ' Compressive sensing predicts that sufficiently sparse vectors can be recovered from 

highly incomplete information using efficient recovery methods such as ^i-minimization. 
Random matrices have become a popular choice for the measurement matrix. In- 
deed, near-optimal uniform recovery results have been shown for such matrices. In 
this note we focus on nonuniform recovery using subgaussian random matrices and £i- 
c/3 I minimization. We provide conditions on the number of samples in terms of the sparsity 

^, ' and the signal length which guarantees that a fixed sparse signal can be recovered with 

a random draw of the matrix using £i-minimization. Our proofs are short and provide 
CNJ ■ explicit and good constants. 

> 

^ ■ 1 Introduction 

(N 

r — I . Compressive sensing allows to reconstruct signals from far fewer measurements than what 

^^ I was considered necessary before. The seminal papers by E. Candes, J. Romberg, T. Tao |5]7j 

and by D. Donoho [TT] have triggered a large research activity in mathematics, engineering 
and computer science with a lot of potential applications. 

In mathematical terms we aim at solving the linear system of equations y = Ax for 
r> I X G C when y G C" and A G (£^xJ^ aj-g given, and when m <^ N. Clearly, in general this 

C^ ' task is impossible since even if A has full rank then there are infinitely many solutions to this 

equation. The situation dramatically changes if x is sparse, that is, ||x||o := i^{£,xi ^ 0} 
is small. We note that || • ||o is called ^o-norm although it is not a norm. 
As a first approach one is led to solve the optimization problem 

min llzllo subject to Az = y, (l-l) 

where Ax = y. Unfortunately, this problem is NP-hard in general, so intractable in practice. 
It has become common to replace the ^o-minimization problem by the ^i -minimization 
problem 

min ||z||i subject to Az = y, (1-2) 

where Ax = y. This problem can be solved by efficient convex optimization techniques [3]. 
As a key result of compressive sensing, under appropriate conditions on A and on the 
sparsity of x, £i-minimization indeed reconstructs the original x. There are basically two 
types of recovery results: 



o 



• Uniform recovery: Such results state that with high probabihty on the draw of the 
random matrix, every sparse vector can be reconstructed under appropriate condi- 
tions. 

• Nonuniform recovery: Such results state that a given sparse vector x can be 
reconstructed with high probability on the the draw of the matrix under appropriate 
conditions. The difference to uniform recovery is that nonuniform recovery does not 
imply that there is a matrix that recovers all x simultaneously. Or in other words, 
the small exceptional set of matrices for which recovery fails may depend on x. 

Uniform recovery via ^i-minimization is for instance satisfied if the by-now classical re- 
stricted isometry property (RIP) holds for A with high probability [HE]. A common choice 
is to take A G ]g"^xj^ as a Gaussian random matrix, that is, the entries of A are independent 
normal distributed mean-zero random variables of variance 1. If 

m>Csln{N/s), (1.3) 

then with probability at least 1 — e~^"^ we have uniform recovery of all s-sparse vectors 
x E M^ using ^i-minimization and A as measurement matrix, see e.g. [71ll3|[T9]. 

In this note we consider nonuniform sparse recovery using Gaussian and more general 
subgaussian random matrices in connection with £i-minimization. Our main results below 
provide non- uniform recovery guarantees with an explicit and good constant. In contrast to 
other works such as J12pi3j we can treat also the recovery of complex vectors. Moreover, we 
get also good constants in the subgaussian case, and in particular, for Bernoulli matrices. 

2 Main results 

2.1 Gaussian case 

We say that an m x N random matrix A is Gaussian if its entries are independent and 
standard normal distributed random variables, that is, having mean zero and variance 1. 
Our nonuniform sparse recovery result for Gaussian matrices and ^i-minimization reads as 
follows. 

Theorem 2.1. Let x € C with ||x||o = s. Let A G ]^™x" ^g ^ randomly drawn Gaussian 
matrix, and let e G (0, 1). // 



m > s 



y/2ln{2N/e) + 1 + \/21n(2/e)/s (2.1 



then with probability at least 1 — e the vector x is the unique solution to the ii-minimization 
problem ()1.2p . 



Remark: In the asymptotic regime A^, s — )• oo, (|2.ip becomes simply 

m> 2sln(3iV/e). (2.2) 



Comparing with (II. 3p we realize that the log-term falls slightly short of the optimal 
one log{N/s). However, we emphasize that our proof is short, and the constant is 
explicit and good. Indeed, when in addition s/N — ;• then we nevertheless reveal 
the conditions found by Donoho and Tanner [121113] . and in particular, the optimal 
constant 2. Note that Donoho and Tanner used methods from random polytopes, 
which are quite different from our proof technique. 



2.2 Subgaussian case 

We generalize our recovery result for matrices with entries that are independent subgaussian 
random variables. A random variable X is called subgaussian if there are constants /3, > 
such that 

F(|^| >t)< /3e~^*' for all t > 0. (2.3) 

It can be shown [22] that X is subgaussian with EX = if and only if there exists a constant 
c (depending only on /3 and 0) such that 

E[exp(AX)] < e^^' for aU A € M. (2.4) 

Important special cases of subgaussian mean-zero random variables are standard Gaussians, 
and Rademacher (Bernoulli) variables, that is, random variables that take the values ±1 
with equal probability. For both of these random variables the constant c = 1/2, see also 
Section 12.31 

A random matrix with entries that are independent mean-zero subgaussian random 
variables with the same constant c in (j2.4p is called a subgaussian random matrix. Note 
that the entries are not required to be identically distributed. 

Theorem 2.2. Let x E C^ with \\x\\q = s. Let A G i^^ixTV ^g ^ random draw of a 
subgaussian matrix with constant c in (12. 4p . and /ei e G (0, 1). // 

m > ^75 sln(4iV/e), (2.5) 

1 - (^) / ln(4iV/e)-i/2in(4/e)V2 

(where we assume additionally that N, e are such that the denominator above is positive) 
then with probability at least 1 — e the vector x is the unique solution to the ii-minimization 
problem ()1.2p . The constant C in (I2.5p only depends on c. 

More precisely, the constant C = 1.646c~^, where c = c{c) is the constant from Lemma 
(jE.ip below, see also Lemma lE. 21 

Remark: If we consider the asymptotic regime A^ — )• oo, the number of measurements that 
guarantees recovery with high probability scales like (4c)sln(A^). 

2.3 Bernoulli case 

We specialize the previous result for subgaussian matrices to Bernoulli (Rademacher) ma- 
trices, that is, random matrices with independent entries taking the value ±1 with equal 
probability. We are then able to give explicit constants for the constants appearing in the 
result of Theorem 12.21 If y is a Bernoulli random variable, then 

E(exp(Ay)) = -{e^ + e"^) < e^^" 

The last inequality can be derived by using Taylor series. This shows that the subgaussian 
constant c = 1/2 in the Bernoulli case. Further, we have the following concentration 
inequality for a matrix B G W^'^^ with entries as independent realizations of ±l/-y/m. 



I RtI|2 _ ||t-||2 
2 2 



>t||x||2)<2e-f(*'/2-*V3)^ (2.6) 



for all X G M^, t G (0, 1), see e.g. [Il[2]. We can simply estimate t^ < t^ in ([Z6]) and get 
c = 1/12 in Lemma IE. 21 and consequently C = 19.76. 



Corollary 2.3. Let x € C with \\x\\q = s. Let A € ]^"^x^ ^g ^ matrix with entries that 
are independent Bernoulli random variables, and let e G (0,1). // 

m > 7 r^-n^ — r^TT7^s\n{4N/e), (2.7) 

- 1-5.45 ln(4iV/e)-V2in(4/e)V2 ^ /^' v ; 

then with probability at least 1 — e the vector x is the unique solution to the ii-niinimization 
problem (jl.2p . 



2.4 Relation to previous work 

Recently, there have been several papers dealing with nonuniform recovery. Most of these 
papers only consider the Gaussian case while our results extend to subgaussian and in 
particular to Bernoulli matrices. 

As already mentioned, Donoho and Tanner [13] obtain nonuniform recovery results 
(terminology is "weak phase transitions") for Gaussian matrices via methods from random 
polytopes. They operate essentially in an asymptotic regime (although some of their results 
apply also for finite values of A'^, m, s). They consider the case that 

m/N —7- 5, s/m — t- p, log(A^)/r7T, — t- 0, A^ — t- cc, 

where /?, 5 are some fixed values. Recovery conditions are then expressed in terms of p and 
6 in this asymptotic regime. In particular, they get a (weak) transition curve pwi^) such 
that p < pwi^) implies recovery with high probability and p > pwi^) rnean failure with 
high probability (as A^ — )• oo). Moreover, they show that pw{^) ~ 21og((5~"^) as 5 — )■ 0. 
Translated back into the quantities N,m,s this gives m > 2slog(A^) in an asymptotic 
regime, which is essentially (|2.2p . 

Candes and Plan give a rather general framework for nonuniform recovery in [8] , which 
applies to measurement matrices with independent rows having bounded entries. In fact, 
they prove a recovery condition for such random matrices of the form m > Cs In(A^) for 
some constant C. However, they do not get explicit and good constants. Dossal et al. [14] . 
derive a recovery condition for Gaussian matrices of the form m > csln(A^), where c 
approaches 2 in an asymptotic regime. These both papers also obtain stability results for 
noisy measurements. 

Finally, Chandrasekaran et al. [9] use convex geometry in order to obtain nonuniform 
recovery results. They develop a rather general framework that applies also to low rank 
recovery and further setups. However, they can only treat Gaussian measurements. They 
approach the recovery problem via Gaussian widths of certain convex sets. In particu- 
lar, they estimate the number of Gaussian measurements needed in order to recover an s 
sparse vector by m > 2s{ln{p/s — 1) + 1) which is essentially the optimal result. It is not 
straightforward to extend their method to subgaussian measurements as they heavily use 
the rotation invariance of Gaussian random vectors. 

3 Proofs 

3.1 Notation 

We start with setting up some notation needed in the proofs. Let [A^] denote the set 
{1,2,..., A^}. The column submatrix of a matrix A consisting of the columns indexed by 



S is written As = {aj)j^s where S C [N] and aj £ M™, j = 1, . . . ,m denote the columns of 
A. Similarly x^ G C'^ denotes the vector x € C^ restricted to the entries in S, and x G C^ 
is called s-sparse if supp(x) = {i : xi j^ 0} = S with S C [A^] and |5| = s, i.e., ||2;||o = s. 
We further need to introduce the sign vector sgn(x) G C^ having entries 

sgn(x), :={ ^ if^i/0, 
^ ^ ^^ I if xj = 0, -^ ^ ^ 

The Moore-Penrose pseudo-inverse of a matrix B such that [B*B) is invertible is given by 
i?t = [B* B)~^ B* , so that i?"''i? = Id, where Id is the identity matrix. 

3.2 The Gaussian case 

We set S := supp(3;), which has a cardinality s. By Corollarv IA.21 for recovery via ii- 
minimization, it is sufficient to show that 

K(^5)V,sgn(x^))| = |(a,,(4rsgn(x^))|<l for all i€[N]\S. 

Therefore, the failure probability for recovery is bounded by 

V := ¥{31 5 K(^5)^a^,sgn(x^))| > 1). 

If we condition X := (a^, (Ag)* sg'n.[x^ )) on ^4^, it is a Gaussian random variable. Further, 
-^ — Z^?Li('^^)j[(^s)*sg^(^'^)]i is centered so its variance u'^ can be estimate by 

m 

u' = nx') = J^E[(a,)2][(4)*sgn(x^)]| 

= ||(4)*sgn(x^)||i < a^l{As)\Hn{x')\\l = a^l{As)s, 

where (Tmin denotes the smallest singular value. The last inequality uses the fact that 
11(4)112^2 = Il4ll2^2 = <^^]ni^s)- Then it follows that 

V<¥(3£^S K(As)^a£,sgn(x^))| > 1 ||(4)*sgn(x^)||2 < a 



+ P(||(4)*sgn(x^)||2>a) 
< 2iVexp(-l/2a2) + P{a;;^l{As)V^ > a). (3.1) 

The inequality in (j3.1|) uses the tail estimate (jC.ip for a gaussian random variable, the 
union bound, and the independence of ai and As- The first term in (jS.ip is bounded by 

e/2if 

a < , (3.2) 

^21n(2iV/e) 

The second term in (j3.1|) can be estimated using (IB.lh below, 

IP«in(^5) V^ > «) 

= P(an,in(^5) < ^sja) = P {araU^s/V^) < V^/iV^a)) 



-m(l - (q^^ + l)^/sJm)'^ 



< exp ^ '-— ^^ ' ' . (3.3) 



If we choose a that makes (j3.2p an equality, plug it into condition (j3.3p . and require that 
(jS.Sp is bounded by e/2 we arrive at the condition 

-l2 



m > s 



V21n(2iV/e) + 1 + V21n(2/e)/s 



which ensures recovery with probability at least 1 — e. This concludes the proof of Theorem 
EH D 



3.3 Subgaussian case 

We follow a similar path as in the proof of Gaussian case. We denote S := supp(x). We 
can bound the failure probability V by 

V<¥(3i^S K(A5)^a^sgn(x^))| > 1 |||(4)*sgn(x^)||2 < a^ 

+ P(||(4rsgn(a;^)||2>a). (3.4) 

The first term in (13. 4p can be bounded by using Lemma ID. II Conditioning on As and 

||(4)*sgn(x'^)||2 < a we get 

m 

F{\{{As)^a,,sgn{x'))\ > 1) = P(| j;(a,),[(4)*sgn(x^)],-| > 1) < 2exp(-l/(4ca2)). 

So by the union bound the first term in (j3.4p can be estimated by 2A^exp(— l/(4ca^)), which 
in turn is no larger than e/2 provided 

a < y^l/{4cln{4N/e)). (3.5) 

For the second term in (j3.4p . we have 

P(||(4)*sgn(x^)||2 > a) < na^l{As)V~s > a) 

= F{amin{As) < V^/a) = P (a^miAs/yM) < ^^ ) . 

Lemma IE. H and Lemma lE?2] implv that a matrix B := Asj ^Jm with normalized subgaussian 
rows satisfy 



P(^mm(5) < \/r^^) < nWB - Id||2^2 > 5) < e/2 

provided m > C5~'^{3s + ln(4e~^)), where C depends on subgaussian constant c. The 

choice —;=— = \/l — S yields 6 = 1 ^. Combining these arguments and choosing a 

that makes (13. Sp an equality, we can bound the failure probability by e provided 

^^^^/ _4c.ln(4iVA)y' 

\ m J 

We define the variable 7 := 1 - {^f^ ln(4iV/e)-^/2 in(4/e)V2. Observe that 7 G (0, 1) for 
A^ large enough. If 7 > 4csln(4A^/e)/7n, that is, if 

?n> — sln(47V/e), (3.7) 

7 

then condition ()3.6p is implied by 

m> 3C7(l-7)-2(s + ln(4/e)/3). (3.8) 

If we plug 7 into (j3.7p and (j3.8p , it can be seen that (j3.7p implies (j3.8p . This completes the 
proof of Theorem 12.21 D 



A Recovery conditions 

In this section we state some theorems that were used in the proof of main theorem, directly 
or indirectly. The proofs of Theorems 12.11 and 12.21 require a condition for sparse recovery, 
which not only depends on the matrix A but also on the sparse vector x G C to be 
recovered. The following theorem is due to J.J. Fuchs |15j in the real- valued case and was 
extended to the complex case by J. Tropp [21], see also [201 Theorem 2.8] for a slightly 
simplified proof. 

Theorem A.l. Let A G £,mxN ^^^ ^ ^ (^N ^^^j^ ^ ._ s^pp(^3;)_ Assume that As is 
injective and that there exists a vector h G C" such that 

A*sh = sgn(x^), 
\{A*h)e\ <1, £e [N]\S. 

Then x is the unique solution to the ii-minimization problem (jl.2p with Ax = y. 

Choosing the vector h = (^^) sgn(x'^) leads to the following corollary. 

Corollary A. 2. Let A G ([^"i.><N ^^^ ^ ^ (^N ^j^^j^ g ._ supp(^j;)_ jf ^/^g matrix As is 
injective and if 

|((^5)'^a^,sgn(x^))| < 1 for all £ G [iV] \ 5 

then the vector x is the unique solution to the ii-minimization problem (|1.1|) with y = Ax. 

B Singular values of Gaussian matrix 

An elegant estimation for the smallest singular value of a normalized Gaussian matrix 
B G M™^'*, where the entries of B are independent and follow the normal distribution 
A/'(0, l/m), was provided in [10], 

P(amin(i?) < 1 - \AM -r)< e-™^'/^ (B.l) 

Its proof relies on the Slepian-Gordon Lemma [16^117] and concentration of measure for 
Lipschitz functions |18] . 



C Tail estimate for a gaussian random variable 

For a mean-zero Gaussian random variable X with variance a^ we have the tail estimate 

Indeed, a mean-zero Gaussian variable g with variance satisfies by [20 1 Lemma 10.2] 

P(|g| > t) = ^ / e-'"/^dt < e-''l\ 
V 2vr Jt 

Rescaling gives the tail estimate ()C.1|) . 



D Tail estimate for sums of subgaussian variables 

The following estimate for sums of subgaussian random variables appears for instance in 
[22]. 

Lemma D.l. Let Xi, . . . , Xm be a sequence of independent mean-zero subgaussian random 
variables with the same parameter c as in ()2.4p . Let a G M.^^ be some vector. Then 
Z := X^fci ^j-^j ^^ subgaussian, that is, for t > 0, 

M 

P(| ^fljXjl >t)< 2exp(-i2/(4c||a|||)). 
Proof. For convenience we provide a proof. By independence we have 

M M M M 

Eexp{e^ajXj) = E JJexp(6'ajXj) = JjEexp((9ajXj) < ]^exp(6'ajXj) 

j=l 1=1 i=l 1=1 

= exp(c||a||2^^). 

This shows that Z subgaussian with parameter c||a||| in (|2.4p . We apply Markov's inequality 
to get 

P(Z >t)= P(exp(eZ) > exp(0t)) < E[exp(^Z)]e~''* < e^il^lla^""^*. 

The optimal choice 6 = i/(2c||a||2) yields 

F{Z>t) < e-*'/(4c|l'^lli). 
Repeating the above computation with —Z instead of Z shows that 

P(-^ > i) < e"*'^^^'"""^\ 
and the union bound yields the desired estimate P(|^| > t) < 2e~* /(^^ll'^lb). □ 

E Concentration Inequalities 

The following concentration inequality for subgaussian random variables appears, for in- 
stance, in [TlfTO]. 

Lemma E.l. Let A be an mx N random matrix with independent, isotropic, and subgaus- 
sian rows with the same parameter c as in (|2.4p . Then, for all x G M^ and every t S (0, 1), 
normalized matrix A = -n=A satisfies 

^{\\\^A\2 - \\A\l\ > A\A\l) < 2exp(-cmt2), (E.l) 

where c depends only on c. 

Combing the above concentration inequality with the net technique we can derive the 
following estimate on the condition of (submatrices of) subgaussian random matrices. While 
this is well-known in principle the right scaling in 5 seemingly has not appeared elsewhere 
in the literature, compare with [21ll9j. 



Lemma E.2. Let S C [N] with card(5) = s. Suppose that m x N random matrix A is 
drawn according to a probability distribution for which the concentration inequality (jE.ip 
holds, that is, for t > 0, 

P(|Px||2 - \\x\\l\ > t\\x\\l) < 2exp(-cmt2) for all x G R^ , 

for some c G IR. Then, for 6 £ (0, 1), 

\\ A*s As -Id\\2^2 <S 

with probability at least 1 — e provided 

m > C6^^{3s + ln(2e-i)), (E.2) 

withC = 1.646c-^ 

Proof. Since most available statements have an additional log((5~^)-term in ()E.2p . we include 
the proof of this lemma for the sake of completeness. 

Let p G (0, v2 — 1) be a number to be determined later. According to a classical covering 
number argument, see e.g. |20l Proposition 10.1], there exists a finite subset U of the unit 
sphere 5 = {j; G M^,supp(j;) C S, \\x\\2 = 1}, which satisfies 

1^1 < I 1 H — ) and minug(/||z — ulU < p for all z £ S. 

V PJ 

The concentration inequality (jE.ip yields 

P(|||^u||2 - \\u\\l\ > t\\u\\2 for some u£K) 
<^r{\\\Au\\l-\\u\\l\ >t\\u\\l) <2\U\exp{-ct 



m) 



ueu 

<2( 1 + - 
P 



j exp (— ct ?7i) . 



The positive number t will be set later depending on 5 and on p. Let us assume for now 
that the realization of the random matrix A yields 

iP^^lli - Il'"ll2| < ^ foralluGC/. (E.3) 

By the above, this occurs with probability exceeding 

1-2M+-J exp{-ct^m). (E.4) 

Next we show that (jE.SP implies |||Ax||2 — Ha^Hil ^ ^ foi" all x G 5, that is HA^A^— Id||2-i.2 < 
6 (when t is determined appropriately). Let B = ^^^4^ — Id, so that we have to show 
||^||2-5.2 < S- Note that (jE.Sp means that \{Bu, u)\ < t for all u £ U. Now consider a vector 
X G 5, for which we choose a vector u £ U satisfying ||x — it||2 < p < \/2 — 1. We obtain 

\{Bx,x)\ = \{B {u + X — u) , u + X — u)\ 

= \{Bu, u) + {B{x — u),x — u) + 2 {Bu,x — u)\ 

< \{Bu,u)\ + \{B (x — u) , X — u)\ + 2 ||i?n||2 ||x — ii||2 

<t+ 11^112^2 P' + 2 11^112^2 P- 



Taking the supremum over all x G 5, we deduce that 



S||2^2 <* + 11^112^2 (P+2p), i.e., \\B\ 



< 



2^2 - 9_ /'^_Ln2 



2-(p + l)^ 



Note that the division by 2 — (p+ 1)^ is justified by the assumption that p < \/2 — 1. Then 
we choose 

t = ts,p := (2 - (p + if) 6, 

so that ||i?||2-5.2 < S, and with our definition of t, 

P {\\A*sAs - Id||2^2 >S)<2(l + -\ exp {-c6\2 - {p + iffm) . (E.5) 

Hence, Hyd^A^ — Id||2-s>2 < ^ with probability at least 1 — e provided 

"^ ^ g(2_J+ 1)2)2 '^'' iH^ + 2/p)^ + ln(2e-^)) . (E.6) 

Now we choose p such that ln(l + 2/p) = 3, that is, p = 2/(e^ — 1). Then ()E.6p gives the 
condition 

m > C5~^ (3s + ln(2e^^)) (E.7) 

with C = 1.646 c~^. This concludes the proof. D 

Acknowledgement 

The authors would like to thank the Hausdorff Center for Mathematics for support, and 
acknowledge funding through the WWTF project SPORTS (MA07-004). We would also 
like to thank Simon Foucart. Indeed, the proof of Lemma lE. 21 is taken from a book draft 
that the second author is currently preparing with him. 

References 

[1] D. Achlioptas. Database-friendly random projections: Johnson-Lindenstrauss with 
binary coins. J. Comput. Syst. Sci., 66(4):671-687, 2003. 

[2] R. G. Baraniuk, M. Davenport, R. A. DeVore, and M. Wakin. A simple proof of the 
restricted isometry property for random matrices. Constr. Approx., 28(3):253-263, 
2008. 

[3] S. Boyd and L. Vandenberghe. Convex Optimization. Cambridge Univ. Press, 2004. 

[4] E. Candes. The restricted isometry property and its implications for compressed sens- 
ing. C. R. Acad. Sci. Paris S'er. I Math., 346:589-592, 2008. 

[5] E. Candes, J. Romberg, and T. Tao. Robust uncertainty principles: exact signal re- 
construction from highly incomplete frequency information. IEEE Trans. Inf. Theory, 
52(2):489-509, 2006. 

[6] E. Candes, J. Romberg, and T. Tao. Stable signal recovery from incomplete and 
inaccurate measurements. Coram. Pure Appl. Math., 59(8): 1207-1223, 2006. 

10 



[7] E. Candes and T. Tao. Near optimal signal recovery from random projections: uni- 
versal encoding strategies? IEEE Trans. Inf. Theory, 52(12):5406~5425, 2006. 

[8] E. J. Candes and Y. Plan. A probabilistic and ripless theory of compressed sensing. 
Preprint, 2010. larXiv:1011.3854[ 

[9] V. Chandrasekaran, B. Recht, P. A. Parrilo, and A. S. Willsky. The convex geometry 
of linear inverse problems. Preprint, 2010. I arXiv: 1012.0621] 

[10] K. Davidson and S. Szarek. Local operator theory, random matrices and Banach 
spaces. In W. B. Johnson and J. Lindenstrauss, editors, Handbook of the geometry of 
Banach spaces I. Elsevier, 2001. 

[11] D. Donoho. Compressed sensing. IEEE Trans. Inform. Theory, 52(4):1289-1306, 2006. 

[12] D. Donoho and J. Tanner. Thresholds for the recovery of sparse solutions via 11 
minimization. In Conf. on Information Sciences and Systems, 2006. 

[13] D. Donoho and J. Tanner. Counting faces of randomly-projected polytopes when the 
projection radically lowers dimension. J. Amer. Math. Soc, 22(l):l-53, 2009. 

[14] C. Dossal, M.-L. Chabanol, G. Peyre, and J. Fadili. Sharp support recovery from noisy 
random measurements by 11 minimization. Preprint, 2011. I arXiv: 1101. 15771 

[15] J. -J. Fuchs. On sparse representations in arbitrary redundant bases. IEEE Trans. Inf. 
Th, page 1344, 2004. 

[16] Y. Gordon. Some inequalities for Gaussian processes and applications. Israel J. Math., 
50(4):265-289, 1985. 

[17] Y. Gordon. Elliptically contoured distributions. Probab. Theory Related Fields, 
76(4):429-438, 1987. 

[18] M. Ledoux. The Concentration of Measure Phenomenon. AMS, 2001. 

[19] S. Mendelson, A. Pajor, and N. Tomczak Jaegermann. Uniform uncertainty principle 
for Bernoulli and subgaussian ensembles. Constr. Approx., 28(3):277-289, 2009. 

[20] H. Rauhut. Compressive Sensing and Structured Random Matrices. In M. Fornasier, 
editor. Theoretical Foundations and Numerical Methods for Sparse Recovery, volume 9 
of Radon Series Comp. Appl. Math. deCruyter, 2010. 

[21] J. Tropp. Recovery of short, complex linear combinations via li minimization. IEEE 
Trans. Inf. Theory, 51(4):1568-1570, 2005. 

[22] R. Vershynin. Introduction to the non-asymptotic analysis of random matrices. In 
Y. Eldar and G. Kutyniok, editors, Compressed Sensing: Theory and Applications. 
Cambridge Univ Press, to appear. 



11 



